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Abstract 

The evolution of a bound system in the expanding background has been investigated in this paper. 
The background is described by a FRW universe with the modified holographic dark energy model, 
whose equation of state parameter changes with time and can cross the phantom boundary. To study 
the evolution of the bound system, an interpolating metric is considered, and on this basis the geodesics 
of a test particle are given. The equation of motion and the effective potential are also derived from the 
geodesics. By studying the the effective potential and the evolution of the radius of a test particle in the 
bound system of the Milky Way galaxy, we have found that the galaxy would go through three stages: 
expands from a singular point; stays in a discoid for a period of time; big rip in the future. With the help 
of analysing the critical angular momentum, we find that the test particle needs less angular momentum 
to escape from the center mass as time passes. 
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1 Introduction 

There is a great deal of observational evidence to prove that our universe is undergoing an accelerating 
expansion phase [3111111 [5]. To explain this phenomenon, a large number of theoretical models have 
been put forward, and among these models the dark energy models have attracted a lot of attention [1 
0 n n [n [n HI- Although numerous dark energy models have been proposed, the nature of it is still 
very largely in the dark. Among these dark energy models, the so called holographic dark energy, which 
arises from the holographic principle [HilllllllS], seems to provide a more simple and reasonable frame 
to investigate the problem of dark energy [H UHl- According to the holographic principle, the energy 
density of the holographic dark energy is written as px = 3c^MpL“^, where c is a numerical constant, and 
Mp = l/CSnG is the reduced Planck mass and L is the largest IR cutoff. Until now, various types of 
IR cutoff have been considered to solve the dark energy problems, such as the Hubble radius [niiniin], 
the particle horizon imiia, the future event horizon [231 [53] , the cosmological conformal time |5S1 [55] , or 
other generalized IR cutoff [53 [551 US] [301 [23 133 1331133] • Nowadays, the researches on the holographic 
dark energy have attracted so many scientists, and lots of remarkable works have been done in this field 

[35ll36ll33l3ai39jl3njl33l4ai43ll4g. 

However, the current accelerating expansion can not be derived from the holographic dark energy models 
with Hubble horizon or particle horizon as the IR cutoff mss]- Meanwhile, Setare et al. 133 SS] proved 
that the holographic dark energy model in the non-flat universe enclosed by the event horizon cannot cross 
the phantom divide line. To solve these problems, Gao, et al. m raised a new holographic Ricci dark energy 
model, where the length scale is given by L ~ and R is the Ricci curvature scalar. Granda, et al. 

[481149] proposed a new form of IR cutoff for the holographic dark energy, which is a linear combination of 
H and H^, i.e., px oc aH^ + pH. This model could avoid the causality problem, which appears when using 
the event horizon area as the IR cutoff, meanwhile, the fine tuning problem and the coincidence problem 
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can also be solved effectively. Recently, Chimento, et al. [50] proposed a modified form of holographic Ricci 
dark energy model, it is 

= ( 1 ) 

where a and /3 are free constants. After their work, the modified holographic Ricci dark energy (MHRDE) is 
studied in connection with the interacting dark matter in many papers, and more information can be found 

in [5T115^ 1551 154] . 

In modern physics, the relation between cosmic expansion and local physics is a problem of principle in 
general relativity which still awaits a definitive answer [53 ISnilSZlISHj- The physics of a bound system in 
the expanding universe is an immediate approach for this problem. A great deal of attention have been put 
on the effects of the universe’s expansion on gravitationally bound system such as planetary, galaxies and 
cluster systems [33 ISHl ISIl EH IMl [Ml ES] • Among them, Nesseris, et al. [5S] have studied different bound 
systems in the phantom and quintessence universe, and they also gave the numerical reconstruction of the 
dissociating bound orbits. In their work, they found that the bound system got unbounded at the time when 
the minimum of the time-dependent effective potential disappeared. This is not the time when the phantom 
energy gravitational potential, which is due to the average (p -I- 3p), balance the attractive gravitational 
potential of the mass M of the system. In the present paper, we would like to examine a bound system in an 
expanding universe filled with dark matter and MHRDE. In such a model, the equation of state parameter 
of MHRDE is time-depending. 

This paper is organized as follows: In Sec. 13 we give a brief review of the MHRDE model and the 
evolution of its equation of state parameter. The evolution of energy densities for dark matter and dark 
energy are studied next. In Sec. [3 using an interpolating metric, which is Schwarzschild like in small 
scales and a general time-depending FRW universe in large scales, we give the test particle’s radial equation 
of motion. Based on this motion equation, the effective potential is given, too. Then, the bound system 
corresponding to the Milky Way galaxy is studied by discussing the evolution of the effective potential and 
the change of the circular orbits. The conclusions are given in Sec. |4| In this work, we assume today’s scale 
factor oq 1; far future a —>■ oo, and the unit SttG = 1. 


2 A brief review of the universe with the MHRDE model 

Now, considering that there are only two components, i.e., the pressureless dark matter and the negative 
pressure dark energy. Then the Friedmann equation describing the evolution of the FRW universe can be 
written as 

= \ {Pm + Px) , (2) 

where pm and px are the energy densities of dark matter and dark energy, H = a/a is the Hubble parameter 
and a is the scale factor. The total energy density must satisfies the following conservation law: 

Pm ff Px T {Pm “t” Px “t” Px) — 0- (3) 

Since we suppose that there is no interaction between dark matter and dark energy, the conservation law 
would be written in the following two separate equations 

Pm + ^Hpm = 0, ( 4 ) 

Px+3H{px+Px) = 0- (5) 

Let’s define that, 

^ ~ TTH' Pm = , Px = , (6) 

where Hq is the present value of the Hubble parameter. Substituting Eq. (jB]) into Eq. ®, we obtain the 
equation for conformal density of dark matter as 

Pm — ^mO^ (’f) 

where = Pmo/^H^ is the present density of dark matter. Using Eqs. (|T]), (jB]) and the Friedmann 
equation can be written in the following conformal form 

= ^moa-^ + 3 (^ + 3ah^) . ( 8 ) 
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The general solution of the above differential equation is given by 

= f5|f|„oa-3 ^ (9) 

where (7 is a constant and can be determined by the initial condition h‘^\x=Q = 1- Using the initial condition 
one can obtain (7=1 — Now, the conformal density of the MHRDE can be identified as 


h = 




,-3/3 


( 10 ) 


and the conformal Hubble parameter is 

a — I, 


= fir. 


-3 


loa 


1 - 


1-/3 




lOfl 


a — p 




mO I ^ 


,-3/3 


Using Eqs. and m, the equation of state parameter for the MHRDE can be written as 


= -l 


(a-l)n^o+/3[l-/3-(a-;3)0,^o]a^(^-'^/ 


( 11 ) 


( 12 ) 


From Eq. (USD, one can find that the equation of state parameter provides the possibility of transition from 
Wx > — 1 to Wa; < —1, which corresponds to the quintom model. If we consider a universe dominated by 
MHRDE, where the contribution from the DM is negligible, then Eq. m becomes Wx = —1 + P- So if 
P < 0, Wx can cross the phantom divide. 

However, from the above equations one can find that there are two free parameters that need to be 
determined. To determine one of them, we will take a special condition into consideration. The derivative of 
H with respect to time is given by 77 = — f where u = Pmipx is the ratio between dark matter 

and MHRDE, then Eq. (fTTl) would be written in the following term 


h^ = 


l-\-U 

ot—0 


ah^ - 


(13) 


From this equation we obtain that the value of /3 given in terms of the free parameter a as /3 = u(l — a) + 
(1 + Wx)- Taking the boundary conditions and the ratio between dark matter and MHRDE, i.e.. 


PmO 

'^xO — 1? ^0 — — Z pr 

PxO J- ^ ^mO 


into consideration, one can obtain 

^ = re(i-«)- (14) 

From Eq. (HU, it is easy to find that /3 decreases as a increases. Now, the free parameters have been reduced 
to one, and for a = 1 we have /3 = 0. 
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Figure 1: The evolution of the state parameter of MHRDE with respect to Ino, we choose flmo = 0.3. 

The evolutions of the state parameter of MHRDE and the densities of dark matter and MHRDE are 
plotted in Figs. [I]and[21 respectively. Combining the two figures, one can easily find that the future of the 
universe would be MHRDE dominated. When a = 1, as shown in Fig. [U the equation of state parameter 
of MHRDE is ma; = — 1, and as shown in Fig. ^ its energy density is a constant, which indicates that the 
MHRDE is ACDM like when a = 1. 
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Figure 2: The evolution of the energy density of dark matter and MHRDE, we choose flmo = 0.3. 


3 Evolution of the bound system 

In order to study the bound system in the expanding universe, one needs the geodesics of a test particle 
in the expanding universe. Now, let’s consider an appropriate metric that satisfies the static Schwarzschild 
metric at small distance and a time-depending FRW space time at large distances. Following most authors 
ElESlEl], we consider such an interpolating metric, which is under the Newtonian limit (weak field, low 
velocities), takeing the form [65] : 


ds^ = 


2GM\ 

{t)p) 


dt^ - aitf ■ [dp^ 

a[t)p/ 

+p‘^ ■ {dO^ + sin^ Odp"^)], 
where p is the comoving radial coordinate. Using r = a{t)p, one can obtain the geodesics 

- - f d - ^ + = 0^ 

= L, 


(15) 

(16) 
(17) 


where L is the constant angular momentum per unit mass. Then, the radial equation of motion for a test 
particle in the Newtonian limit would be 


r = — 


GM 


+ -r. 


(18) 


From this equation one can obtain the effective potential, which determines the dynamics of the bound 
system, by f = —dVjdr. Integrating it with respect to r and ignoring the arbitrary integration function of 
time, one can obtain the effective potential 


Here, taking ^ into consideration, we can obtain 


r = ^ = Hr' 

where r' = dr/dlno, then for r one has 

We also take d/a = H + into account and use Eq. m, then, Eq. (IlSp would be 

GM 


(19) 

( 20 ) 

( 21 ) 


W'-fi(i7^)V= + ^ 

Similarly, one can obtain the effective potential in terms of H as 

TA _ GM I G 1 rl /■ E72'>/ I rr2l„2 

- - ) +H \r . 


( 22 ) 


(23) 
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In what follows, we will focus on the evolution of a two body bound system in the expanding unverse. 
Let us assume that at the present time, i.e., a = oq, a test particle moves on a circular orbit with radius ro, 
and — GM/tq. Taking the following conformal forms into consideration 

Wo - r - r/ro, 14// - Veff/{H^rl), 

(24) 

Eqs. (1331) and (1331) could be written in the following forms 


/iV'+ i(h2)V =-^ + ^ + A(a)2r, 

(25) 

Fe// = -^ + ||-iA(a)V, 

(26) 

where 


A(a)^ = + h^. 

(27) 

By analysing the effective potential, one could find the radius of the circular orbit is 
T'min (o) of the effective potential. From Eq. (1261) , we know that the location of rmin 
scale factor. Radius of the circular orbit is given by solving the following equation 

1 given by the minimum 
(a) is depending on the 

„ 1 1 _ n 

' min ' mzn -r ± u. 

(28) 


In fact, Eq. (1^ has a solution only for A(a)^ < [5S]. Therefore, when A(a)^ = the 

minimum of the effective potential disappears, then, the system becomes unbound. Here, we define that the 
solution of equation A(a)^ = 27a;Q/256 as a = Oi-ip, which means the big rip occurs at the trip time Orip. 

To investigate the future radial evolution of a bound gravitational system, we will take a specific bound 
system — the Milky Way galasy (M = 2 x ro = 5 x lO^^cm, wq = 182) — into consideration. The 

numerical evolution of the radius for different a is shown in Fig. [31 As shown in Fig. O it is easy to find that 
all the curves have two turning point, one occurs at small a, which means the early universe, and another 
occurs at big a, which means the future universe. According to Fig. |3l we find that r < tq occurs before the 
first turning point, and r > vq occurs after the second turning point. From this we can say that the Milky 
Way galaxy would go through three stages: the first stage is the expanding stage: from a very small radius 
to today’s size; the second is the steady stage: the radius changes little and the galaxy stays in a discoid; the 
third is the rip stage: the radius of the galaxy becomes bigger and bigger that the galaxy would be ripped. 
From Fig. |3l one can also find that as a increases the values of ari-p decreases, which means that the big rip 
occurs more early. The two body bound system corresponding to the Milky Way galaxy is plotted in Fig. 
m From Fig. 21 it is more easy to find that the galaxy comes from a singular point, stays in discoid for a 
period of time, and will be ripped finally. 



ln(a/ao) 


Figure 3: The numerical evolution of the radius for the Milky Way galaxy is plotted under different cases: 
continuous line for a = 1.1, dashed line for a = 1.2 and dotted line for a = 4/3. Here we choose 12^0 = 0.3. 

The evolution of the effective potential around the trip time arip is shown in Fig. |5l As a contrast, the 
effective potential of the present time is also plotted in Fig. [SI As we know, the stability of a bound system 
depends on whether the effective potential has a minimum. From Fig. [SI one can find that when a < arip, 
as the radius r increases the effective potential curves decreases to a minimum first, then increases to a 
maximum and finally decreases to —oo, and its maximum decreases as the scale factor increases. Then we 
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Figure 4: The evolution of the system of the Milky Way galaxy until the potential minimum disappears. 
Here, we choose a = 1.2 and Hmo = 0.3. 


know that the effective potential curves have a minimum, which indicates the possibility of bound orbits for 
the test particle. However, as we can see from Fig. [5l when a = o^ip, the effective potential curve has only 
one inflection point and when a > Urip there is no inflection point. Then, we can say that when the scale 
factor a > Grip, the gravitational force can not bound the test particle any longer, they will fly to distant 
area under the negative pressure force of the MHRDE. 



Figure 5: The evolution of the effective potential corresponding to different scale factor a with respect to 
the dimensionless radius r/rp, and we choose a = 1.2, Umo = 0.3. 


We have studied the test particle at a circular orbit above, where we assume that = GM/tq. 

Now, let’s consider a different case that the test particle has different angular momentums. Eq. ([^ could 
be written in the following form 


K// = -^ + ^-iA(a)2 


(29) 


where = L’^KH^Tq). From Eq. (l29ll one can find that when = |a;o[wo/(4A(a)^)]^, the effective 

potential curves have no extremum. The critical angular momentum Lc with respect to In a is plotted in Fig. 
m From Fig. ini one can obtain that the critical angular momentum decreases as In a increases, which means 
with the passage of time, the test particle would need less angular momentum to escape the gravitational 
force from the center mass. 

In order to make this more clear, the effective potential curves for the test particle with different angular 
momentums are shown in Fig. [71 From Fig. [7] (a), one can find that the critical angular momentum L is 
more big than ujq, i.e., >> Wg. Fig. [7] (b) shows that = Wq. As we all know, when the effective 

potential has no extremum, the test particle could not be bounded together. Then, with the help of the 
two graphics in Fig. [71 we can say that the test particle would need less angular momentum to escape the 
immense pull from the center mass with the passing of time. 
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Figure 6: The critical angular momentum with respect to In a. Here, we choose a = 1.2 and = 0.3. 




(a) 


Figure 7: The evolution of the effective potential with respect to In a under different angular momentum. 
In figure (a), we choose the scale factor as a = oq and in figure (b) we choose a = Urip- Here, we choose 
a = 1.2 and Hmo = 0.3. 


4 Conclusion 

In this paper, we have investigated the evolution of a bound system in the expanding universe, the universe 
is described by the FRW universe with only two cosmic components: dark matter and dark energy. In the 
present paper the dark energy is considered as the MHRDE, whose density is px = and the 

equation of state parameter of it is time-depending. From Sec. [2] one can find that the MHRDE behaves 
like the quintom. 

To investigate the bound system in the expanding universe, we have considered an interpolating metric 
which can reduce to the static Schwarzschild metric at small scales but ERW metric at large scales. The 
geodesics of a test particle are derived by this metric and using the geodesics we get the radial equation of 
motion for a two body bound system. The test particle’s motion in the Milky Way galaxy is examined by 
analysing the equation of motion and the effective potential. We have found that the the galaxy would go 
through three stages: the first stage is expanding from a singular point, then it will stay in a discoid for a 
period of time, and it will be ripped owing to the negative pressure of the MHRDE hnally. 

Moreover, by studying the evolution of the critical angular momentum with respect to the scale factor 
a, we find that with the passage of time, the test particle would need less angular momentum to escape the 
immense pull of the center mass. 


Acknowledgments 

This work was supported in part by the National Natural Science Foundation of China (Grants No. II347I0I, 
No. II405076), the Science Foundation of the Education Department of Yunnan Province (Grant No. 
20I4Y066), and the Talent Cultivation Eoundation of Kunming University of Science and Technology (Grants 
No. KKSY20I207053, No. KKSY20I356060). Yu Zhang would like to acknowledge the support of the 
working Eunds of the Introduced High-level Talents of Yunnan Province from the Department of Human 
Resources and Social Security of Yunnan Province. 


7 























References 

[1] Riess, A. G., et al.: Astron. J. 116, 1009 (1998) 

[2] Perlmutter, S., et al.: Astrophys. J. 517, 565 (1999) 

[3] Knop, R. A., et al.: Astrophys.J. 598, 102 (2003) 

[4] Tonry, J. L., et al.: Astrophys. J. 594, 1 (2003) 

[5] Riess, A. G., et ah: Astrophys. J. 607, 665 (2004) 

[6] Sahni, V. and Starobinsky, A. A.: Int. J. Mod. Phys. D 09, 373 (2000) 

[7] Padmanabhan, T.: Phys. Rep. 380, 235 (2003) 

[8] Peebles, P. J. E. and Ratra, B.: Rev. Mod. Phys. 75, 559 (2003) 

[9] Sahni, V. and Starobinsky, A.: Int. J. Mod. Phys. D 15, 2105 (2006) 

[10] Gopeland, E. J., Sami, M. and Tsujikawa, S.: Int. J. Mod. Phys. D 15, 1753 (2006) 

[11] Frieman, J., Turner, M. and Huterer, D.: Annu. Rev. Astron. Astrophys. 46, 385 (2008) 

[12] Li, M., Li, X. D., Wang, S. and Wang, Y.: Commun. Theor. Phys. 56, 525 (2011) 

[13] ’t Hooft, G.: arXiv:gr-qc/9310026 (1993) 

[14] Susskind, L.: J. Math. Phys. 36, 6377 (1995) 

[15] Fischler, W. and Susskind, L.: arXiv:hep-th/9806039 (1998) 

[16] Bousso, R.: Rev. Mod. Phys. 74, 825C874 (2002) 

[17] Setare, M. R.: Phys. Lett. B 642, 1-4 (2006) 

[18] Sheykhi, A.: Phys. Lett. B 680, 2 (2009) 

[19] Gohen, A. G., Kaplan, D. B. and Nelson, A. E.: Phys. Rev. Lett. 82, 4971 (1999) 

[20] Horava, P., Minic, D.: Phys. Rev. Lett. 85, 1610 (2000) 

[21] Thomas, S.: Phys. Rev. Lett. 89, 081301 (2002) 

[22] Bousso, R.: J. High Energy Phys. 07, 004 (1999) 

[23] Huang, Q. G. and Li, M.: JGAP 0408, 013 (2004) 

[24] Li, M.: Phys. Lett. B 603, 1 (2004) 

[25] Gai, R. G.: Phys. Lett. B 657, 228 (2007) 

[26] Wei, H. and Cai, R. G.: Phys. Lett. B 660, 113 (2008) 

[27] Ghen, S. and Jing, J.: Phys. Lett. B 679, 144 (2009) 

[28] Gai, R. G., Hu, B. and Zhang, Y.: Gommun. Theor. Phys. 51, 954 (2009) 

[29] Duran, 1. and Pavon, D.: Phys. Rev. D 83, 023504 (2011) 

[30] Wang, Y. and Xu, L.: Phys. Rev. D 81, 083523 (2010) 

[31] Xu, L., Li, W. and Lu, J.: Eur. Phys. J. C 60, 135-140 (2009) 

[32] Xu, L., Lu, J. and Li, W.: Eur. Phys. J. C 64, 89-95 (2009) 

[33] Xu, L. and Wang, Y.: JGAP 6, 002 (2010) 


8 




[34] Zhang, M. J., Ma, C., Zhang, Z. S., Zhai, Z. X. and Zhang, T. J.: Phys. Rev. D 88, 063534 (2013) 

[35] Bouhmadi-Lopez, M. and Tavakoli, Y.: Phys. Rev. D 87, 023515 (2013) 

[36] Cai, Y. F., Saridakis, E. N., Setare, M. R. and Xia, J. Q.: Phys. Rep. 493, 1C60 (2010) 

[37] del Campo, S., Fabris, J. C., Herrera, R. and Zimdahl, W.: Phys. Rev. D 87, 123002 (2013) 

[38] Ghaffari, S., Dehghani, M. H. and Sheykhi, A.: Phys. Rev. D 89, 123009 (2014) 

[39] Li, M., Li, X. D., Meng, J. and Zhang, Z.: Phys. Rev. D 88, 023503 (2013) 

[40] Setare, M. R. and Vagenas, E. C.: Int. J. Mod. Phys. D 18, 147 (2009) 

[41] Setare, M. R.: Phys. Lett. B 644, 99-103 (2007) 

[42] Setare, M. R.: Eur. Phys. J. C 50, 991-998 (2007) 

[43] Setare, M. R.: Int. J. Mod. Phys. D 17, 2219 (2008) 

[44] Xu, L.: Phys. Rev. D 87, 043525 (2013) 

[45] Hsu, S. D. H.: Phys. Lett. B 594, 13 (2004) 

[46] Setare, M. R.: Phys. Lett. B 642, 421-425 (2006) 

[47] Gao, G., Wu, F., Chen, X. and Shen, Y. G.: Phys. Rev. D 79, 043511 (2009) 

[48] Granda, L. N. and Oliveros, A.: Phys. Lett. B 669, 275 (2008) 

[49] Granda, L. N. and Oliveros, A.: Phys. Lett. B 671, 199-202 (2009) 

[50] Chimento, L. P. and Richarte, M. G.: Phys. Rev. D 84, 123507 (2011) 

[51] Chattopadhyay, S. and Pasqua, A.: Indian J. Phys. 87, 1053-1057 (2013) 

[52] Chimento, L. P. and Richarte, M. G.: Phys. Rev. D 85, 127301 (2012) 

[53] Chimento, L. P., Forte, M. and Richarte, M. G.: Eur. Phys. J. C 73, 2285 (2013) 

[54] Chimento, L. P., Forte, M. and Richarte, M. G.: Eur. Phys. J. C 73, 2352 (2013) 

[55] Bars, L, Steinhardt, P. and Turok, N.: Phys. Rev. D 89, 043515 (2014) 

[56] Faraoni, V. and Jacques, A.: Phys. Rev. D 76, 063510 (2007) 

[57] Gao, C., Chen, X., Faraoni, V. and Shen, Y. G.: Phys. Rev. D 78, 024008 (2008) 

[58] Gao, C., Chen, X., Shen, Y. G. and Faraoni, V.: Phys. Rev. D 84, 104047 (2011) 

[59] Bona, C. and Stela, J.: Phys. Rev. D 36, 2915 (1987) 

[60] Bonnor, W. B.: Mon. Not. R. Astron. Soc. 282, 1467 (1996) 

[61] Cooperstock, F. L, Faraoni, V. and Vollick, D. N.: Astrophys. J. 503, 61 (1998) 

[62] Einstein, A. and Straus, E. G.: Rev. Mod. Phys. 18, 148 (1946) 

[63] Gonzalez-Diaz, P. F. and Siguenza, C. L.: Phys. Lett. B 589, 78 (2004) 

[64] Stefancic, H.: Phys. Lett. B 595, 9 (2004) 

[65] Nesseris, S. and Perivolaropoulos, L.: Phys. Rev. D 70, 123529 (2004) 

[66] McVittie, G. C.: Mon. Not. R. Astron. Soc. 93, 325 (1933) 


9 



